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Abstract 

Context. The high Galactic latitude sky at millimeter and submm wavelengths contains significant cosmological information about 
the early Universe (in terms of the cosmic microwave background) but also the process of structure formation in the Universe from 
the far infrared background produced by early galaxies and the Sunyaev-Zeldovich effect in clusters of galaxies. 
Aims. As the Planck mission will produce full sky maps in this frequency range, deeper maps of selected low-foregrounds patches of 
the sky can produce complementary and important information. Here we analyze the performance of a balloon-borne survey covering 
a 10° x 10° patch of the sky with a few arcminute resolution and very high pixel sensitivity. 

Methods. We simulate the different components of the mm/submm sky (i.e., CMB anisotropics, SZ effect, radio and infrared sources, 
far infrared background, and interstellar dust) using current knowledge about each of them. We then combine them, adding detector 
noise, to produce detailed simulated observations in four observational bands ranging from 130 to 500 GHz. Finally, we analyze the 
simulated maps and estimate the performance of the instrument in extracting the relevant information about each of the components. 
Results. We find that the CMB angular power spectrum is accurately recovered up to I ~ 3000. Using the Sunyaev-Zel'dovich effect, 
most of the galaxy clusters present in our input map are detected (60% efficiency overall). Our results also show that much stronger 
constrains can be placed on far infrared background models. 

Key words. Cosmic microwave background - Methods: data analysis - Infrared: galaxies - Submillimeter 



1. Introduction 

The high Galactic latitude sky at far infrared and millimetric fre- 
quencies is mainly dominated by four components: the cosmic 
microwave background (CMB), the Sunyaev-Zel'dovich (SZ) 
effect in clusters of galaxies, Galactic dust, and the far infrared 
background (FIRB) from unresolved point sources, as well as 
radio and infrared sources. The "cosmological window" extends 
roughly from 70 to 200 GHz, consisting at lower frequencies 
of the interstellar emission of spinning dust grains, free-free and 
synchrotron emission from the interstellar medium, which domi- 
nate over the cosmological microwave background, and at higher 
frequencies of the clumpy foreground from "cirrus clouds" of 
interstellar dust and the background signal due to unresolved 
sources at high redshift, which dominate over the sky brightness. 

The CMB provides, together with its anisotropies, signifi- 
cant information on the primordial phases and the evolution of 
the Universe. Its frequency spectru m can be distorted by dif- 
ferent effects ( Agh anim et al.ll2008h . occurring on both local or 
global scales: in this work, we only consider the local distor- 
tion due to the Sunyaev-Zel'dovich effect caused by the interac- 
tion between the hot gas present in the clusters of galaxies and 
the photons of the cosmic background radiation. To allow the 
extraction of the cosmological parameters, we need to separate 
these d ifferent types of emission, hence to use multi-band experi- 
ments dJuin et al.l2007l) . Moreover, the study of these foreground 
signals is very useful to improve our knowledge about struc- 



ture formation of the primordial Univ erse dGispert et alJ f2000: 
lAghanim et d]|2005tlPires et al.ll2006l) . 

In this paper, we develop an accurate simulation of the ex- 
pected signals in the 1 30-500 GHz frequency range, and demon- 
strate how it is possible to separate and recover the input sig- 
nals. We use data from a multi-band, high resolution survey of 
10° X 10° of the sky, which is feasible using a long-duration 
balloon telescope. As baseline, we assume the characteristics of 
the OLIMPO experiment. The paper is organized as follows: in 
Sect. [2] we briefly summarize the characteristics of OLIMPO; in 
Sect.|3]we describe the astrophysical components included in our 
model; Sect. [4] explains how we modeled each component and 
the map-making process; Sect. [5] describes the map-separation 
algorithm we used; finally, we discuss our results in Sect. [6]and 
report our conclusions in Sect. [7] 



2. The OLIMPO experiment 



Send offprint requests to: luca.conversi@esa.int 



OLIMPO is a balloon borne experiment (Masietal ] 120051) 
whose purpose is to observe selected portions of the sky on four 
frequency bands centered on v = 143, 217, 353, and 450 GHz, 
with an angular resolution ranging from 5' FWHM at 143 GHz 
to 2' FWHM at 450 GHz (see Table [TJ>. It is constituted by a 
Cassegrain telescope of 2.6 m of diameter, advanced bolomet- 
ric detectors, and a long duration 3 He cryostat, which will be 
mounted on a stratospheric balloon payload able to perform long 
duration (10 - 15 days) circumpolar flights. 

The incoming radiation is focused onto four focal planes 
using a system of three dichroics: each focal plane consists of 
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an array of TES (transition edge superconductor) detectors. The 
four arrays are composed of 19 detectors in the two lower fre- 
quency bands (143 and 217 GHz) and 24 in the others. Each ar- 
ray fills the optically correct area of the focal plane (~ 15' in di- 
ameter projected onto the sky). The cryogenic reimaging optics 
is mounted inside the experiment's cryostat, cooled to 2 K, while 
the bolometers are cooled to 300 mK: extensive baffling and 
a cold Lyot-stop reduce both straylight and sidelobes by about 
30 dB. The optical system was optimized to allow diffraction- 
limited operations even with significant tilt of the primary mirror 
during sky scans. Radiation is refocused into the cryogenic op- 
tics through a side of the cryostat, to allow operation of the sys- 
tem with the telescope axis horizontal or at low elevation, during 
ground-based calibrations or in-flight observations of planets. 

The high angular resolution and the wide range of frequen- 
cies covered by OLIMPO allow us to achieve the following sci- 
entific goals: 

- measurement of the CMB anisotropies at high multipoles 
(( > 2000): by taking advantage of its wide frequency cov- 
erage and the high sensitivity, OLIMPO will significantly 
improve our current knowledge, allowing the study of the 
damping tail of the angular power spectrum of the CMB; 

- measurement of the Sunyaev-Zel'dovich effect caused by 
cluster of galaxies: it is not possible to study the charac- 
teristics of high redshift clusters using X-ray data, since the 
signal is too faint to be detected in a blind survey with the 
present X-ray satellites. The Sunyaev-Zel'dovich effect, in- 
stead, is independent of the redshift of the cluster, thus a fun- 
damental tool for cosmological measurements. The survey of 
known clusters and the deep blind observation of a blank sky 
area with OLIMPO will permit a more complete statistical 
analysis of the clusters distribution and of their evolution; 

- study of the far infrared background (FIRB) due to unre- 
solved galaxies at high redshift: this measurement will im- 
prove our knowledge of the galaxy formation process, and in 
general of the epoch immediately following recombination. 

Because of its larger telescope, OLIMPO has a higher angu- 
lar re solution than Planck HFI in the same bands dLamarre et alj 
l200l . Moreover, working as an observatory, it can integrate 
longer than Planck on selected sky regions. This, combined to 
the larger number of detectors per band, results in a higher sen- 
sitivity for the selected targets. 

OLIMPO also co mplements B LAST dPascale et alJ 120081) 
and Herschel SPIRE (Grif fin et all 1201 Ol) . covering lower fre- 
quency bands and thus allowing clearer component separation, 
especi ally for the SZ effect; on the other hand, i t complements 
ACT dFowlerll2004 and SPT dRuhl etalj|2004l) at higher fre- 
quencies and with comparable angular resolution. 



Channel 


143 


217 


353 


450 


Central frequency [GHz] 


146 


221 


360 


458 


Bandwidth [GHz FWHM] 


40 


61 


98 


126 


Field of view radius [arcmin] 


20 


20 


20 


20 


Resolution [arcmin FWHM] 


5.2 


3.7 


2.3 


1.9 


Number of detectors 


19 


19 


24 


24 


Detector's NEP \jaK/ Vflz] 


137.2 


240.7 


316.4 


533.5 



Table 1. Characteristics of the OLIMPO instrument in the 4 
wavebands. Filters' properties, detector's NEP, and resolution 
come from laboratory measurement of flight samples. Each de- 
tector array fills the same unvignetted, circular field of view. 



3. Model parameters 

As described in Sect. [1] the sky brightness in the frequency 
range 130 - 500 GHz is mainly dominated by four compo- 
nents: cosmic microwave background, Sunyaev-Zel'dovich ef- 
fect, Galactic dust, and far infrared background. In addition to 
these signals, both instrumental noise and the contribution of ra- 
dio and infrared point sources must be counted. Hence, we built 
a model to describe the expected signals in four frequency bands 
(nominally centered on v = 143, 217, 353, and 450 GHz) to 
check the capabilities of the OLIMPO experiment, and whether 
it can separate the different effects and to what extent. 

3.1. Cosmic microwave background 

The c osmic microwave background (CMB, Penzia s & Wilson! 
Il965l) is the brightest background of the millimeter sky. This 
radiation and its properties are the strongest co nfirmations of the 
standard cosmological model dGamowlll948l) . It has a thermal 
black -body spectrum at a temperature of 2.725 K dMather et al.l 
1999), thus the spectrum peaks at 160 GHz, while the anisotropy 
spectrum peaks at 220 GHz. At present, the more robust mea- 
surement available of the angular distribu tion of the CMB i s 
that derived from the WMAP 5 -year data (Dun klev etal1 l2009): 
using its estimation of the CMB angular power spectrum, it is 
possible to generate CMB simulated maps as Gaussian random 
realizations of the Fourier transform of the CMB spectrum. 

3.2. Sunyaev-Zel'dovich effect 

The SZ effect is the result of the interaction between the CMB 
photons and the hot gas of electrons (T e > 5 keV) in clus- 
ters of galaxies by inverse Compton scattering, i.e. the photons 
adsorb part of the electron's ener gy: this leads to a blueshift 
of the CMB freq uency spectrum (Zeldovich & Sunvaev 1969; 
iBirkinshawll 19991) . Two contributions can be identified: a ther- 
mal one, directly related to T e , and a kinetical one, which is re- 
lated to the cluster peculiar velocity V r . 

In the non-r elativistic limit, al so known as the Kompaneets 
approximation dKompaneetsll 19571) . the frequency spectrum of 
the thermal SZ effect is described by 
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where x — hvlk^T cmb and the comptonization parameter y 
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where the integration is performed along the line of sight. The 
thermal SZ effect is described by the product of the analytic 
expression of the spectrum and the comptonization parameter. 
Thus, in this approximation, the spectrum does not depend on 
the cluster's parameters. 

The comptonizatio n parameter is obtained from hydrody- 
namical simulations (Borga nTet al.l l2005t ISpringel et al. 2005) 
that follow the process of formation and evolution of galaxy 
clusters in a cosmological framework. These simulations include 
an advanced treatment of several processes, such as radiative 
cooling, star formation, thermal conduction, release of metals, 
and energy feedback from supernovae type la and II. The re- 
sult of these simulations is a box containing a very high number 
of pa rticles (~ 10 10 in the largest one, the Millennium simula- 
tion, Springel 2005): by integrating along the line of sight, it is 
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possible to obtain an accurate simulation of a map of the comp- 
tonization parameter. 

As already mentioned, when a cluster has a non-zero pecu- 
liar velocity in the CMB reference frame, there is a kinematic 
effect in addition to the purely thermal SZ one. In the reference 
frame of the scattering gas, CMB radiation appears anisotropic, 
and the effect of the inverse Compton scattering is to slightly 
re-isotropize the radiation. 

Assuming that the thermal and kinematic effects are inde- 
pendent, it is possible to calculate the kinematic correction to 
the spectrum 



(Il998h . SCU BA dDve et 



AI K (x) 



2 (kg Tcmb) 
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where T e — cr T J n e dl is the optical depth. In principle, a pho- 
tometric measurement at x - 3.83 (v 220 GHz) is able to 
detect a non-zero cluster's peculiar velocity along the line of 
sight: given this parameter, it is possible to subtract the kine- 
matic SZ signal at the other frequencies and distinguish the ther- 
mal component. Therefore, this effect can be used to measure 
the large-scale motion of objects at high redshifts. However, the 
spectrum of the kinematic SZ is identical to the spectrum of pri- 
mary CMB anisotropy, so it is almost impossible to distinguish 
the two contributions along the same line of sight. Therefore, in 
the following, we do not consider this effect. 

3.3. Dust 



The Schlegel, Finkbeiner and Davis model #8 (Schlege l et al.l 
1 1998b is an accurate description of the far infrared and mi- 
crowave emission from the diffuse interstellar du st in the Galaxy- 
It is b ased on high resolution, far infrared IRAS (Wheelock et al, 
11994 data and calibrated with the DIRBE (Bogg ess et al.lll992l) 
data at v = 3000 GHz, which have a much higher calibration 
accuracy but poorer spatial re solution. T his model then ext ends 
to FIRAS dFixsen etal.ll 1994 and DMR (ISmoot et al.ll 19901) fre- 
quencies, 100-2100 GHz and 3 1 .5 - 90 GHz respectively. The 
outputs of the model are full sky maps at submillimeter and mil- 
limeter wavelengths. 

This is a multicomponent model, because it takes into ac- 
count an interstellar medium of many types of molecules and 
dust grains, and fits all available data in the range 200 - 
3000 GHz. The model analyzes only the sky regions where the 
far infrared emission is expected to be dominated by diffuse in- 
terstellar medium. Hence, a spatial mask that excludes all the 
undesirable regions is created, leaving a map covering ~ 7 1 % of 
the sky observed by FIRAS. 

This model was developed for the limit of quite big grains, 
0.01 < a < 0.25 fim, where a is the grain radius: they reach a 
steady thermal equilibrium with the surrounding radiation field. 
Furthermore, since the grain's emission wavelength scales pro- 
portionally to its volume a 3 , larger grains dominate the submil- 
limeter emission. Smaller grains, a < 0.01 /mi, can be heated 
temporarily to very high temperatures and dominate the emis- 
sion at A < 100 fim (v > 3000 GHz), but they neither contribute 
significantly to submillimeter emission, nor are important in the 
FIRAS frequency interval. 

3.4. Far infrared background 

The cosmic infra red backgro u nd (C IB) was first detected in 
COBE map s by [P uget et al. (1996) and l ater endorsed by 
iFixsen et a l. ( 1996). IFixsen et all (I19981 and lBurigana~& Popa 



and BLAST dDevlin et 



tall 



2006), Spitzer (lLagache et al.l2007l) . 
2009) also confirmed its detection. 



The CIB is a product of the cosmological structure formation 
process: matter aggregates into stars and galaxies, then the evo- 
lution of these systems is related to gravitational and nuclear 
processes, which lead to the emission of energy. Cosmic ex- 
pansion and absorption of radiation em itted by dust transform 
this radiative en ergy into CIB emission dHauser & Dw ek 2001; 
iBlain et al.ll200l . 

One of the main problems of far infrared observations is the 
difficulty in resolving sources (confusion noise): this is due to the 
low spatial resolution of the instruments, caused by diffraction 
in small telescopes. Thus, the superimposition of signals com- 
ing from faint unresolved sources form the so-calle d far infrared 
background. Measurements obtained by BLAST dDevlin et al.l 
2009) have been less affected by this problem because this tele- 
scope was able to resolve most of the sources of the CIB. 

Different models have been developed to reproduce 
the FIRB measurements ( e.g., iFranceschinT et al.1 11994 



Pearson & Rowan-Robinson 1996; Burigana etalJ 119971: 
Franceschini et al.ll 1997HGuiderdoni et al.lll998tlXu et al.lr2001 
IValiante et al.l 12009). We adopted the model o f Lagache, Dole 
& Puget (hereafter LDP, lLagache et al.l 120031) . which has the 
advantage of being simple and with a low number of parameters 
and components, but is cabable of reproducing all the available 
observations. In a simple way, it constrains the galaxy luminos- 
ity function evolution with redshift, and fits the existing source 
number counts and redshift distributions, and the CIB intensity 
and its fluctuations, from the mid to the far infrared wavebands 
dFernandez-Conde et al . 2008). This model is based on both 
local spectral templates of starburst and normal galaxies and the 
local infrared luminosity function at 60 fim. 

The LDP model distinguishes the FIRB sources into two cat- 
egories, starburst and normal galaxies, and empirically differen- 
tiates between their evolutions. AGNs dominate the emitted ra- 
diation only at very high luminosities (L > 2 ■ 10 12 L G ) and have 
a negligible influence at the wavelength of interest for OLIMPO. 

Given a source at redshift z, its flux S v at a given wavelength 
A — Aq can be written as a function of the rest frame luminosity 
L v (expressed in W/Hz) 



S v (L,z, A = A Q ) = 



(l+z)K(L,z)L v (A = Ao) 



AnD\ 



(4) 



where Di is the luminosity distance and K(L,z) is the K- 
correction factor, defined as 



K(L,z) = 



L v (z=0) 

The number of sources per solid angle and redshift interval is 



(5) 



-^- ( L,z) = M,(L,z) ( l + z)^, 
dz a log L dz 



(6) 



where dV/dz is the differential element of volume, once the cos- 
mological parameters are fixed, N is the number of sources per 
unit volume and luminosity interval as a function of redshift, and 
No is given by the luminosity function. 

Differential and integrated counts for a given flux S and at 
wavelength A-Aq can be written as 

dN dz 

(L,z)— (L,z)d log L, (7) 



-= f f 

dS J L J : 



z dzdlogL dS 
C dN 



(8) 
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Hence, the FIRB intensity Ifirb produced by sources with flux 
S < Smo.x, expressed in Jy/sr, is 



'FIRB 



(*S Max 

Jo 



S ds dS > 



(9) 



while the fluctuation densities caused by the sources below the 
observational flux limit S o, which correspond to both the con- 
fusion noise and sensitivity limit of the experiment, are given 
by 

Pfluc= I ' S '— < 10 > 



Jo S dS' 



thus Pfluc is expressed in Jy 2 /sr. In the case of OLIMPO, the 5- 
<t confusion noise So ranges from 1.4 to 11.5 mJy in the four 
wavebands. 



3.5. Point sources 

The contribution of point sources to the total signal is mainly that 
of radio sources in the two lower frequency bands, and resolved 
infrared sources (the unresolved ones being already taken into 
account in the FIRB model) in the higher frequency bands. Their 
contribution to the total signal is negligible with respect to the 
CMB at low frequencies and the FIRB at higher frequencies. In 
particular, the r.m.s. of the fluctuations due to radio sources is 4% 
of the r.m.s. due to the CMB in the OLIMPO beam at 143 GHz 
and 3% of the r.m.s. due to the FIRB in the OLIMPO beam at 
450 GHz. However, it is important to include these sources in 
the simulation, because their presence may alter the results of 
the component separation (see Sect. [3j and can strongly affect 
the estimate of the angular power spectra (see Sect. 16.31 ) at high 
multipoles, with an increase of up to 20% at i ~ 2000. 

We modeled these two kind of sources with a differential 
source count dN(S)/dS of the form 



dN(S) 
dS 



A N s-y< 

A N S 



72-71 
Brk ' 



for S > S Br k, 

~ 72 for S < S Brk- 



(11) 



For the radio sources, we used the WMAP catalogue 
dWright et al.l [2009) to estimate the parameters of Eq. [TTJ at 
143 GHz, obtaining yi =0.23,y 2 = 1.03,5^ = 1.9xlO- 3 mJy, 
and the normalization factor An = 59. 

For the infrared sources, we fitted the dN(S)/dS from 
the LDP model with the previous analytical expression 
(iFriedmann & Bouc het 2004), deriving y x = 3.3, y 2 = L8, 
S Brk = 1.9 mJy, and A N = 18 • 10 9 . 



4. Map generation 

We generated 10° x 10° maps, by taking into account the in- 
strument characteristics and using as a baseline the OLIMPO 
experiment (see Sect. [2}. A beam of 5' FWHM was used for 
all the wavebands. This means that we degraded the map angu- 
lar resolution for the high frequency channels, so that all maps 
are immediately comparable. We assumed arrays of 19 detectors 
in the two lower frequency bands (143 and 217 GHz) and 24 
detectors in the others. We decided to integrate for 100 h over 
the observed area of the sky, which is centered on the Galactic 
coordinates b = 65°, / = 155°. These assumptions are con- 
sistent with a flight of ~ 15 days departing from the Svalbard 
islan ds (Norway) in June -July and circumnavigating the north 
pole dPeterzen et al.ll2008l) . 



4.1. Cosmic microwave background map 

The CMB map is a Gaussian random realization with a null av- 
erage value and an angular power spectrum that was measured 
by WMAP and is defined by the cosm ological parameters of the 
ACDM model dDunklev et alj|2009t) . We generated a full sky 
map using the synfa st program (part of the HEALPix package of 
iGorski et al.l (|2005)): it was generated for nside - 2048, which 
roughly corresponds to a pixel dimension of 5/3'. Hence, each 
map's pixel is identified by the Galactic coordinates (latitude b 
and longitude €) and its flux is Acmb- Then, using the smooth- 
ing program included in the HEALPix package, the map is con- 
volved with the instrument beam working in harmonic space. In 
Galactic coordinates (b, €), the convolved map S cmb is given by 

S C MB(b,0= f f ' A CMB (b',(')f (7) cos b'db'df, (12) 

Jo J-jt/2 

where AcMB(b', £') is the input CMB map towards the Galactic 
coordinates (b' , £'), 



/(r) 



i 



2ncri 



exp 



H 


( 7 )1 




\ O" FWHM j 



(13) 



FWHM 

is a symmetric Gaussian angular response with o- FWHM = 
5'/ V8ln2 and 

y - arccos(sinb sinfo' + cosb cosb' cos({ - (')) (14) 

is the angular distance between direction (b, I) and direction 
(b',C). Finally, we extracted a 10° x 10° portion of sky. The 
final map used for the analysis is presented in the top-left panel 
of Figure [TJ 

4.2. Sunyaev-Zel'dovich effect map 



We calculated the SZ spectrum in t he non-relativistic limit (i.e. 
in the Kompaneets approximation, Kompa neetsl 1 1 9571) . In this 
case, the spectrum of all the clusters present in the map is in- 
sensitive to the gas temperature: as shown by Eqs. [TJand[2l the 
comptonization parameter is indeed frequency-independent and 
provides only spatial information about the clusters. 

Hence, to simulate the SZ effect, we produced a map of the 
comptonization parameter y(b, () using hydrodynamical simula- 
tions. The SZ spectrum A$z was then computed in the four ex- 
periment's wavebands and converted into K 

2(k B T CMB ) 4 f 00 IP 6 T C MBE B (v)g(v) 



dB(v,T) 
dT 



■ dv 



= 10 6 T, 



CMB 
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Jo 



E B (v) ■ 



hv 



kT, 



■ coth 



CMB 



hv 



2kT t 



CMB 



T-Tcmb 



dv. 



where 



g(x) = 



(e 



-') ? |.v.^i-4| 



with x - 



hv 



kT, 



(15) 



(16) 



CMB 



and E B (v) represents the pass-band filter spectrum of the four 
channels. We then multiplied A$z by the map of the comptoniza- 
tion parameter y(b, €). Finally, the map is convolved with the 
beam as done for the anisotropy map (see Eq.lTZb 

SszAb,r>= f [' A sz , B y(b'J')f(y)cosb'db'df. (17) 

J() J-n/2 

Figure [TJ top-right panel, shows the SZ map obtained at 
143 GHz. 
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Figure 1. Input maps of the four components at 143 GHz. Signals are expressed in fiK. At this frequency, apart from the largest 
clusters that produce an intense Sunyaev-Zel'dovic effect, the main contribution is due to CMB anisotropics. 



4.3. Galactic dust map 

As described in Sect. |4] we selected a sky region centered on 
the Galactic coordinates b = 65°, I = 155°: this area is rich in 
known clusters, observable during the whole flight and the dust 
contribution is expected to be very low. 

Dust emission is estimated using the Schlegel, Finkbeiner & 
Davis model #8 (see Sect. I3.31 >. which is able to rescale the IRAS 
100 fim map to the desired frequency. It only requires as input 
the Galactic coordinates and the frequency, giving as output the 
sky temperature in fjK: Ao U st(b, t, v). It is then convolved with 
the filter transfer function and the instrument beam f{y) (see 
Eqs.[l3]and[Tl) 

r*2n /"vr/2 

S Dust , B {b,€) = f(y) cos b'db'df ■ 

Jo J-n/2 

Jr»oo 
A Dus ,(b',t',v)E B (v)dv. (18) 
o 

Figure Q] bottom-left panel, presents the Galactic dust map at 
143 GHz. 

4.4. FIRB map 

Our FIRB maps were compiled in two stages. In the first step, 
we produced four maps in the four OLIMPO wavebands using 



a Gaussian distribution, N(A.FiRB,o~FiRB,b,P), with the average 
and the standard deviation provided by the LDP model in Jy/sr 
and Jy/ respectively. The total signal was then converted 
into /iK. 

In the second step, we randomly distributed the radio and 
infrared point sources over the maps, following the differen- 
tial source counts dN{S)/dS described in Sect. 13.51 Each ra- 
dio source was scaled to higher frequency using a power law 
f(v) oc v~ g , where the spec tral index a varies between -2.5 
and 0.5 (Sokasian et al. I I200T1) . For the infrared sources, we dis- 
tributed only sources with a flux 3-ctfirb above the FIRB back- 
ground, to consider only sources that are not confused (hence 
already taken into account as part of the background). 

Finally, each map is convolved with the instrument beam: 
the resulting map at 143 GHz is shown in Figure [T] bottom-right 
panel. We note that we did not include any clustering component 
in our simulations: this certa inly affects the input FIRB signal 
(iGonzalez-Nuevo et al.l 12005). but not the results discussed in 
Sect. 16.31 Adding the clustering would increase the FIRB power 
spectrum at high multipoles, but the ability to disentangle the 
different components remains the same (and so the ability in re- 
covering the input spectrum), as long as a pixel-based map ex- 
traction method is used (see Sect. [5]l. 
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Figure 2. Total maps predicted at the four frequency bands of OLIMPO: these are the sum of the CMB anisotropies, SZ effect, dust, 
and FIRB components plus the radio and infrared sources and the expected noise in the four bands. Signals are expressed in //K. 
We note that the CMB anisotropy component dominates the 143 and 217 GHz channels, while the 353 and 450 GHz channels are 
dominated by Galactic dust. The far infrared background contribution is comparable to the Galactic dust one, but seeing it is more 
difficult since it lacks of structures at scales larger than beam size. However, when adding the point sources, these loom over the 
other components, particularly at 450 GHz. 



4.5. Noise map 

In building our model, we made no attempt to simulate system- 
atic effects, and we took into account only white noise that has a 
null average and standard deviation equal to 

NET bN pi xe l 

A N oi se ,B = j (19) 

■\jT[ntNBolo,B 

where NETg is the bolometer's noise in fiK/ ^JHz and ranges 
from ~ 150 yuK/ at 143 GHz to ~ 550 fiK/ yftfz at 450 GHz, 
NboIo.b is the number of receivers in each channel, Tj„, is the 
integration time, and Npi xe i is the number of pixel in a row (or 
column) of the map. Each map's pixel has a dimension equal to 
FWHM/3, thus 

N Pi «, = 3 ■ 60' ■ deg = 3^_ 

5' • V3600s- deg 2 5 yT 

where we have taken into account that we are integrating for 
1 h over a 1 deg 2 region. As a consequence, Eq. [19] could be 
rewritten as 

_3 NET B 

&Noise,B - T , — ■ 

3 -\JNboIo.B 



(20) 



The computed noise is finally convolved with the instrument 
beam as done in Eqs.[T2land[T7l 



5. Map extraction 

The total map for each frequency band is obtained by summing 
the maps of each component (see Figure [TJ: the results are the 
four maps shown in Figure [2] 

The signal caused by Sunyaev-Zel'dovich effect is always 
subdominant: the cosmic microwave background anisotropies 
dominate the two lower channels, while at 353 and 450 GHz 
the far infrared background and the Galactic dust contribu- 
tions are of comparable strength. FIRB is not evident in the 
maps because it lacks structure on scales larger than the beam 
size. However, when adding individual point sources, these 
stand out with respect to of the other components, particu- 
larly at 450 GHz. A comprehensive comparison o f the differ- 
ent co mponent-separation methods is reported in iLeach etaD 
However, we decided to use a simple x 2 approach to en- 
sure a tighter control on each step of our analysis. This method 
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Figure 3. Extracted maps of the four components at 143 GHz. Signals are expressed in fiK. The anisotropy map is the most robustly 
extracted; the SZ map clearly shows the brighter clusters; the dust and FIRB component cannot be easily distinguished because 
their spectra are too similar. 



only requires an a priori knowledge of the component frequency 
spectrum but not of their spatial distribution. 

The total signal can be expressed in terms of the input signal 
and its frequency spectrum using the matrix formalism 

S T (&,0 = M-S c (b,€) + Sti(b,(), (22) 

where St is a vector containing the total maps in the four fre- 
quency bands, M is a 4 x 4 matrix whose columns comprise the 
scaling factors for each of the four components relative to one 
frequency band of choice, Sc is a vector consisting of the four 
component maps computed at the chosen frequency, and Sn con- 
tains the four noise maps. 

To disentangle the four components and obtain their maps 
at the chosen frequency, we used the x 2 method. The optimal 
estimation of the input s ignals was obtained b y minimizing, for 
each pixel, the quantity ( iDupac & Giard 200^) 

X 2 = (S T - M ■ Sc) 1 " ■ N-'(S T - M ■ S c ), (23) 

where N is the noise covariance matrix, which in our case is 
reduced to a diagonal matrix whose elements are the noise vari- 
ances at the four frequencies. The solution to Eq.l23lis 

S c (b, () = [M^Mp 1 M^IST 1 S T (b, €), (24) 

where Sc is our optimal estimation of Sc, i.e. it contains the 
signal of the component extracted from the four total maps. 



6. Results 

While the CMB and SZ frequency spectra are theoretically and 
experimentally well constrained, this is not true for the other 
components. As a consequence, we decided to extract the CMB 
and SZ signals and leave the dust, FIRB, and sources in a resid- 
ual map. The results of the component separation are shown in 
Figure[3] where we chose to extract the components at 143 GHz. 

6.1. CMB 

We computed the angular power spectrum on the ex tracted CMB 
map u sing anafast, part of the HEALPix package (Gorsk iet al.l 
2005). However, appropriate corrections mu st be made becaus e 
we observed only a small fraction of the sky dHivon et al. l2002h : 
these include those for instrumental beam shape, scanning strat- 
egy, and instrumental noise. Since some of these characteristics 
will be known only once the experiment have flown, we used 
a simpler approach. We computed the power spectra of both the 
input map and a full sky CMB one, using the same parameters to 
generate them: the ratio of the two power spectra then provided 
us with the correction factor for all the spectra computed in the 
partial maps. 

Figure|4]shows the angular power spectrum used to generate 
the input CMB map (dot-dashed line) compared to the spectrum 
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computed on the extracted CMB map (triangles), using a bin- 
ning of At = Anjbuca = 72, where b/^ax represents the map 
size. Thanks to OLIMPO's characteristics (high resolution, low 
noise receivers, and wide frequency range), we are able to fully 
recover the spectrum up to I 2200; at higher multipoles, our 
data still fit the input spectrum, but the beam contribution is more 
prominent increasing the error in the C( estimation. 
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Figure 4. CMB angular power spectrum used to generate the in- 
put map (dot-dashed line) compared to the recovered spectrum 
(triangles): we used a binning of AC = An /bMax = 72, bMax being 
the map size. 



We also compared the number of clusters found per binned 
flux in both maps, reporting our results in Figure [6] the binning 
is twice the previous standard deviation, i.e. AT = 15.8 fiK. We 
recovered all the clusters brighter than ~ 160^/K, while we lost 
about 50% of the fainter ones due to the lower signal-to-noise 
ratio. 



so r 




50 100 150 200 250 500 

Temperature (^K) 



Figure 6. Comparison of the number of clusters found per 
binned flux in the input (solid line) and extracted (dot-dashed 
line) SZ maps. The binning is twice the standard deviation of 
the flux difference distribution, i.e. AT = 15.8 fiK. 



6.2. SZ effect 

We identified galaxy clusters in the input and extracted S Z maps, 
using the SExtractor algorithm dBertin & Arnoutsll 19961) . choos- 
ing a threshold a-jhr = 3 cr sz , where cr sz is the estimated stan- 
dard deviation of the SZ map's background. We found ~ 270 
clusters in the input map and ~ 170 in the extracted one, without 
any false detection (i.e. clusters found in the extracted map but 
not present in the original one). We then computed the difference 
between the fluxes of the clusters found in both maps: these data 
have a distribution similar to that of a Gaussian centered on zero 
and with a standard deviation <x = 7.9 /jK (see Figure|5j. 



40 7 
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Figure 5. Flux difference of the cluster found in the input and 
in the extracted SZ maps. Data dispersion is well fitted by a 
Gaussian distribution centered on zero and with standard devia- 
tion cr = 7.9 /iK. 



6.3. Dust & FIRB 

The extracted map contains the contribution of Galactic 
dust, FIRB, and point sources. Point sources must be re- 
moved otherwise the angular power spectra can become dis- 
torted at high multipole s. We achieved this using SExtractor 
dBertin & Arnouts| | l996) with a threshold cr T i, r = 3 o- FIRB , where 
o~firb is the standard deviation of the extracted map's back- 
ground as estimated by the program itself. We then computed 
the angular power spectrum of the input dust and FIRB maps, 
separately, and of the extracted map: the left plot of Figure Q 
shows the dust spectrum (yellow line), the FIRB spectrum (black 
line), and the spectrum obtained for the extracted map (blue di- 
amonds). The plot shows that dust is dominant on large scales 
{( < 200), but becomes negligible at higher multipoles. 

We then performed the same analysis after changing the 
FIRB model. In this way, we tested whether OLIMPO is able 
to discriminate between the different evolutionary models: 

- we considered one redshift interval Az — 1 in the range z 

■ 7 and increased the number of sources with respect to the 
original LDP model: this simulated a higher star formation 
rate in this redshift interval; 

- we considered one redshift interval Az = 1 in the range z 

7 and increased the luminosity of sources with respect to 
the original analysis, while the number of sources remained 
the same. 

In the first case, we modified the number of sources, their lu- 
minosities, and their fluxes: increasing the number of sources in 
a given redshift interval indeed caused the total luminosity and 
flux contribution from the same interval also to be augmented. 

We tried different increments, focusing our tests in the range 
Z = 4 — 5, where the star formation rate is not well known. For 
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Figure 7. Angular power spectra obtained on the input maps of dust and FIRB and on the extracted maps. On the left, the yellow 
line is the dust spectrum computed on the input dust-only map; the black line represents the FIRB spectrum computed on the input 
map, using the original LDP model. The blue diamonds are the binned spectra computed on the extracted map. On the right, the 
comparison between the original LDP model and the modified ones when incrementing by 30% the luminosity or the number of 
sources in the redshift range z = 4 — 5. The dashed line illustrates the effect of increasing the luminosity of sources; the green line 
corresponds to an augmentation of 30% in the number of sources in the same redshift range; the solid line is the original FIRB 
spectrum. The yellow and red diamonds are the extracted spectrum for the two modified models, while the blue diamonds are for 
the original one. The error bars represent the 1-cr uncertainty in the computed binned spectra. 
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Figure 8. Comparison between the analytic form of th e infrared 
background measured by FIRAS dFixsen et al.l 11998b and the 
modified FIRB models produced by incrementing by 30% the 
luminosity (crosses) or the number of sources (boxes) in the 
redshift range z = 4-5. The solid and dashed lines represent 
the most accurate estimation and the error at 1-cr in the FIRB 
spectrum obtained with FIRAS; the four boxes and crosses are 
the fluxes of the modified FIRB model at 143, 217, 353, and 
450 GHz. 

each increment in luminosity or number of sources, the corre- 
sponding modified FIRB signal was computed again in each fre- 
quency band, added to the other components, and extracted to- 
gether with the Galactic dust. We computed once more the an- 
gular power spectra of both the modified FIRB input map and 
the extracted one. We then estimated our ability to disentangle 
the different models calculating the Bayes factor Korig/Mod - 
p(.C{\Mo r i S )l p(C(\M M od), where p(C(\M x ) is the marginal like- 
lihood for model M x (i.e. the original one Morig or a modified 
one Mmoci) and is given by 



p{C ( \M x ) oc exp 



I 



(c*-c 



Exp \ 



(25) 



The proportional sign is indicated because of the absence of the 
normalization factor, which vanishes when computing the ratio, 
C^ xp is the binned angular power spectrum obtained on the ex- 
tracted map, while CJ is the one computed on the input FIRB- 
only map. 

The smallest increment for which we were able to disentan- 
gle the original and modified FIRB models corresponds to an 
increase of 30% in the luminosity and of 20% in the number 
of sources. When comparing C^ xp (using as input the original 

model) to C°" g and Cf od , we obtained Kong/\3L = 35 when aug- 
menting the luminosity and Korig/x.iN = 87 when incrementing 
the number of sources. Since the Bayes factors is higher than 10, 
there is strong support for the original LDP model, which was 
indeed the input spectrum of C, p . 

The right plot in Figure|7]shows the results obtained after in- 
creasing by 30% both the luminosities (dashed line and orange 
diamonds represent original and extracted spectra respectively) 
and the number of sources (green line and red diamonds, original 
and extracted spectra respectively) in the range z = 4 — 5, com- 
pared to the original model (solid line and blue diamonds, orig- 
inal and extracted spectra, respectively). In all cases, we were 
able to recover the correct input spectrum. 

We note that the modified models are consi stent with the lim- 
its imposed by FIRAS data ( Fixsen et al . 1998) in all the redshift 
intervals at a 1-cr level: Figure [8] reports the case of increasing 
by 30% both the luminosity (crosses) and the number of sources 
(squares) in the range z = 4 — 5. 



7. Conclusions 

We have presented accurate simulations of the expected 
sky brightness in the OLIMPO wavebands. Using a simple 
component-separation method, we have efficiently disentangled 
the input signals. 

The CMB angular power spectrum is efficiently recovered 
up to 6 - 3000: even if OLIMPO is not specifically designed to 
study the CMB anisotropies, this gives us an idea of its capabil- 
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ities. Moreover, the achieved result allows us to place stronger 
constrains on the damping tail of the CMB spectrum. 

Using the Sunyaev-Zel'dovich effect, we have been able to 
extract most of the galaxy clusters present in our input map: we 
are complete down to 160 fiK and detect 50% of the fainter clus- 
ters, giving a 60% efficiency overall. 

The separation of the dust and FIRB components was ineffi- 
cient at the pixel level. However, after removing the radio and in- 
frared point sources, the computation of the angular power spec- 
tra on the extracted maps was found to be in agreement with the 
sum of the two input components. Moreover, the dust's spectrum 
quickly declines dramatically and its contribution becomes neg- 
ligible on scales smaller than ~ 1° (€ > 200). Even if we had not 
included any clustering component, which should increase the 
input FIRB spectrum at high multipoles, this would not affect 
our ability to disentangle the different components and recover 
their input spectra. 

We were also able to place stronger constrains on far infrared 
background models: we were able to distinguish between the 
original LDP model and a modified one, in which we had in- 
creased the number of sources by 20% or their luminosities by 
30% in a given redshift interval Az = 1 . 

Hence, thanks to its characteristics (high angular resolution, 
several detectors each having low noise, long duration flight cov- 
ering a selected patch of the sky with a long integration time), 
OLIMPO will be able to achieve its scientific goals, comple- 
menting on the one hand the shallower survey of Planck, and on 
the other the lower frequency surveys of bolometric arrays on 
large-size ground-based telescopes. 
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